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Abstract 



This paper discusses the enumeration of two-terminal series-parahel net- 
works, i.e. the number of electrical networks built with n identical elements 
connected in series or parallel with two-terminal nodes. They frequently oc- 
cur in applied probability theory as a model for real networks. The number 
of networks grows asymptotically like R^/n", as for some models of sta- 
tistical physics like self-avoiding walks, lattice animals, meanders, etc. By 
using a exact recurrence relation, the entropy is numerically estimated at 
R = 3.560 839 309 538 943 3 (1), and we show that the sub-leading "univer- 
sal" exponent a is 3/2. 
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I. INTRODUCTION 



The series-parallel networks are one of the simplest models of real networks. Then even 
if it is at last one century old (Mac Mahon 1892), it still occurs frequently in both theoretical 
or applied studies. Some questions like the reliability (Feo and Johnson 1990), the survival 
probability after failures, the routing problem and the queuing delay (Tamir 1993), etc. 
are now an important branch of probability theory applied to communication, computer or 
power networks. 

As usual in this kind of model, the elementary components (resistances or switches) are 
chosen all identical. Then the interest is in the behavior of the structure made of series or 
parallel connections. 

Some authors (Riordan 1958, Lomnicki 1972) investigated this problem with a combina- 
torial point of view, and, in particular, give recurrences to compute the number of different 
structures built with n elements. Then, they consider the problem as formally "solved" . 

In this article, we study this model by using analogies with statistical physics models like 
self-avoiding walks (Des Cloizeaux and Jannink 1987), lattice animals and percolation (Der- 
rida and De Seze 1982), folding and meanders model (Di Francesco et al. 1996), etc. In all 
these models, the number of configurations for n elements (steps, links, folds) grows asymp- 
totically like R^/rf where R depends on the details of the model and a is an "universal" 
exponent, common to a broad set of models ("universality" class). 

In Section 2, we define the problem. In Section 3, we explain how to enumerate the 
series-parallel networks. Finally, in Section 4, we give by numerical extrapolation a value 
for R and we show that a = 3/2. 

II. DEFINITIONS AND NOTATIONS 

The two-terminal series-parallel networks are electric network considered geometrically. 
They are built with elements Xi which are conductor (Xj = 1) or insulator (Xj = 0). When 
two structures (with conductivities Xa and X^) are in parallel, their (global) conductivity 
is Xa + Xb where -|- is the boolean operator Or. When they are in series, it is X1X2 where 
the product is the boolean operator And. 

A network with n elements, which connects two-terminal nodes, has a boolean represen- 
tation F = F{Xi, . . . ,X„). We consider that two networks are equivalent if their boolean 
representations differ only by a permutation of variables {Xj}. For example, XiX2(X3-|-X4), 
Xi(X2 -|- X3)X4 and (Xi -|- X2)X3X4 are equivalent. We note that this definition is less re- 
strictive than topological definition. In other words, equivalence consists of interchange 
inside a group of structures in series, or in parallel. The networks for n = 1, 2, 3, 4 are given 
in Table |. 

Let us call a„ the number of different networks with n elements. It is also the number 
of non-equivalent boolean functions of n variables, built only with And and Or. 

Except for the trivial case n = 1 with oi = 1, we see in Table | two classes of networks : 
essentially series (a-networks) and essentially parallel (n-networks) . There is a duality : by 
changing the boolean And into Or and vice versa, each a-network is dual of a vr-network. 
Let us call 6„ the number of cr-networks (or of vr- networks) . We have 
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TABLE I. Two-terminal series-parallel networks for n = 1.2,3,4. On each row, the two net- 
works are dual. The a„'s (number of networks) are 1, 2, 4, 10, . . .The 6^'s (number of essentially 
series or, by duality, essentially parallel networks) arc 1, 1, 2, 5, ... 



number of essentially 
elements series 




essentially 
parallel 




number of 
networks 








Ai 


1 


2 • • • 


X1X2 




X1+X2 


2 


3 • • • • 


X 1X2X2 

Xi{X2+X3) 




Xi + X2 + A'3 

Xi + A2 A3 


4 




00 



XiX2{X3+X^) 

Xi{X2 + Xs + X^) 

Ai(X2+A3A4) 

{Xi+X2){Xs + Xi) 




X1+X2 + X3 + X4 

Xi + X2 + A'3X4 

Xi + X 2X3X4 

Xi+X2(A3+A4) 

A^iA^2 + X3X4 
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On — 26„ — 5n,l- 



(1) 



III. ENUMERATION AND RECURRENCE GENERATING FUNCTION 



In this section, we explain how to enumerate the different networks, i.e. compute the 
Qn or bn- By use of the duahty, we consider only a cr-network Sn with n > 2 elements. It 
is composed with k 7r-structures of w 1, W2 • • • , Wfc elements, with k > 2 and J2i f^i — n^- Its 
boolean representation is Sn = Pw^Pw^ . . . Pw^. By grouping together the equal Wj, we can 
use a symbolic exponential notation Sn = P^^Pp . . . P^'^. 

The {rj} are non- negative (r^ > 0) integers ; they form a partition of n with at least two 
parts : n — J2j j fj is partitioned in A; = parts. We note a partition as V^2'^'^ . . .q'^". 

For example, the partitions of 5 are {1^, 1^2, 1^3, 12^, 14, 23}. 

Then, to obtain 6„ (the number of o"-network wdth n elements), we first must compute 
the contribution of a given partition, and finally, we will sum over all the partitions. Thanks 
to the above decomposition, the number of Sn with a given partition {rj} is a product of as 
many factors as j, dependent only on bj and rj. 

For example, it is obvious that, for 65, the contribution of the partition (23) is 62^3; 
because 62 is the number of choices for the first 7r-structure of two elements, and 63 for the 
second of three elements. 

When two 7r-structures have the same number of elements, i.e. rj = 2, their contribution 
is not b^ but, to remove duplications obtained by permutations, bj{bj + l)/2. More generally, 
for arbitrary rj > 0, the factor is 
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bj + r J - 1 



T ■ / 

^ / l<Ul<U2...'^Ur - <bj 

J ■' 



E 1- (2) 



So, the bn (or the a„) can be computed by recurrence 



bn = T.I[(^'^"' '1 (3) 



where the sum runs over all the partition {rj} of n with at least two parts, as explained 
above. For example, 65 = 1 + 62 + &3 + ^2(^*2 + l)/2 + &4 + &2&3 = 12. 

With a computer, the bn can be calculated by recurrence. Unfortunately, the number 



of partitions of n grows like exp{ny2n/3)/n (Hardy-Ramanujan formula). We will now 
describe how to obtain a more practical recurrence formula. 
We will use the generating functions 

= 5]6,f = t + t' + 2t3 + 5t^ + ..., (4) 

i>l 

a{t) = ait' = t + 2t^ + At^ + lOt^ + . . . = 2b{t) - t. (5) 

i>l 

Following an argument similar from that used by Euler for the generating function for the 
partition of n, by using Eq. (H) 



^1 \ ^ , Ml + M ^2 , hi + 2 



X . . . 
X . . . 

= l + Ef^;V" + EM" = l + a(t). (7) 

n>l \ / n>2 

As a generic line of the l.h.s. of Eq. (H) is equal to (1 — t-')^*^ , we obtain a formula attributed 
to Mac Mahon (1892) : 

l + a{t) = \[{l-t^r'^ (8) 

i>i 

Eq (1) can be rewritten as a non-local functional relation 

H^ + a{t)] = Y.T.\b^i'' = T.W)- (9) 

k>i i>i k>i 

We obtain a useful recurrence taking the derivative of Eq. (0) 



TABLE II. The a„ for n < 45 



n 




a-n 


n 






an 


n 














an 


1 




1 


16 


4 


507 


352 


31 






306 


481 


637 


180 


676 


2 




2 


17 


14 


611 


576 


32 




1 


039 


825 


579 


149 


516 


3 




4 


18 


47 


633 


486 


33 




3 


533 


238 


105 


431 


144 


4 




10 


19 


156 


047 


204 


34 




12 


022 


695 


376 


246 


364 


5 




24 


20 


513 


477 


502 


35 




40 


964 


814 


174 


562 


504 


6 




66 


21 


1 696 


305 


728 


36 




139 


754 


932 


965 


689 


486 


7 




180 


22 


5 623 


993 


944 


37 




477 


353 


256 


699 


920 


544 


8 




522 


23 


18 706 


733 


128 


38 


1 


632 


304 


839 


666 


178 


522 


9 


1 


532 


24 


62 408 


176 


762 


39 


5 


587 


602 


984 


719 


422 


092 


10 


4 


624 


25 


208 769 


240 


140 


40 


19 


146 


480 


336 


548 


243 


722 


11 


14 


136 


26 


700 129 


713 


630 


41 


65 


670 


379 


223 


051 


730 


824 


12 


43 


930 


27 


2 353 386 


723 


912 


42 


225 


448 


266 


880 


936 


115 


440 


13 


137 


908 


28 


7 927 504 


004 


640 


43 


774 


644 


024 


342 


635 


183 


112 


14 


437 


502 


29 


26 757 247 


573 


360 


44 


2 663 


899 


039 


964 


678 


479 


766 


15 


1 399 


068 


30 


90 479 177 


302 


242 


45 


9 168 


056 


380 


423 


173 


365 


752 



^ a;(t) ^ (10) 

Defining = where the sum runs over the divisors of i including j = 1 and j = i, 

2 naX = (1 + E «.^') E + ^i) (11) 

n>l j>l i>l 

Defining e[ = = — ica, we can compute the a„ with the recurrence formula 

(Lomnicki 1972) 

n-l 

na„ = 1 + + E «n-i(l + ^i)- (12) 
1=1 

for which the computational effort grows like rP. By programming this recurrence on a 
computer with the use of real*16 (i.e. real numbers described by 128 bits), it is possible 
to obtain several thousands of terms with 30 decimal digits of precision. It is more than 
enough to determine the asymptotic behavior. 

This integer sequence belongs to the encyclopedia of Sloane (1996). The beginning is 
given in Table 



IV. ASYMPTOTIC BEHAVIOR 

A brief investigation of Table |l| shows that a„+i ^ 3.5 a„. Then, we can define a 
non- vanishing entropy per element In an/n ~ InR. More generally, it is expected that the 
exponential dominant behavior is corrected by sub-leading power-law for large n 
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TDn 

an - C- (13) 

where C is a constant. 

In fact, in Eq. (p!2D, for large n, e„ ~ na„, then a„ ~ J2i=i ^i^n-i- This last relation 
is identical apart from a shift (n ^ n + 2) to the definition c„+i = J2]=o ^j^n-j of the 
Catalan numbers c„ = {2n)\/[n\{n + 1)!]. Using the Stirling formula, we see that c„ have 
the asymptotic behavior of Eq. (|13D with i? = 4 and a = 3/2. The fact that the a„ and c„ 
have the same asymptotic relation do not imply that they have the same R. As a > 0, in 
the sum, the extremal terms are dominant : coc„ > ciCn-i > C2C„_2 > . . .. Then the value 
of R strongly depends of the value of the first c„ (n small) and the exponential increase of 
ttn and c„ are not described by the same R. On the other hand, we expected the same a as 
we see below. 

This behavior is common in statistical physics. The enumeration of random walks, self- 
avoiding walk, animals, etc. of n step (or links) drawn on a regular lattice always exhibits 
a exponential increase i?" where R depend on the connectivity of the lattice, and a sub- 
leading ?T,~", where a is "universal", i.e. is dependent on the physical nature of the problem 
(dimension, phases, . . . ) but independent on the details of the lattice. 

Let us take a big two-terminal series-parallel network. At each scale, it is "decorated" 
by small parts 6„ (ra = 1, 2, . . .). Then, the contribution of these small 6„ is a large fraction 
of the total entropy InR ; they are similar to the connectivity of the lattice and govern 
the dominant behavior i?". Nevertheless, the sub-leading n~°' depends on more general 
properties and a can be see as an "universal" exponent. It described the singularity of the 
generating function. 

We suppose that the a„ have the following large n expansion 

a„~C:^(l + ^ + ^ + ...) (14) 

and we want estimate the values of R and a by numerical extrapolation. It is more convenient 
to study Vn and Wn 

Vn^ln^ = lnR-^ + ^ + ^ + ... (15) 

U/ _ 2i anan-2 Wi W2 , . 

Wn = In — 5 = a-\ \ 7 + . . . (16) 

By plotting Vn (resp. Wn) versus 1/n, the hypothesis (|T^ is consistent if the data draw a 
straight line ; its intersection with the y-axis gives Ini? (resp. a). 

A more systematic scheme of extrapolation is the Richardson polynomial extrapolation 
process (Brezinski and Redivo Zaglia 1991). Supposing Eq. (|14|), 

^AP(nPVn) =lnR + 0(1/^^+^) (17) 
p] 

^AP(nPWn) =a + 0(l/nP+i) 
p] 

where A is the difference operator A/„ = /„ — This method gives the limits (i? and a) 
and show whether the hypothesis (0) is consistent or not. When p is large, the acceleration 
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FIG. 1. The generating function a{t). Its singularity occurs at t* = 1/R w 0.28. At this point, 
a{t*) = 1 and the function has a square root behavior because a = 3/2. 

of convergence is faster ; unfortunately, the numerical instabilities are amplified and after a 
given n, extrapolated values (Eq. |1^) become chaotic. In this work, better extrapolations, 
obtained with p between 3 and 5, are 

R = 3.560 839 309 538 943 3 ± 10"^^ (18) 
a = 1.5±10"^l (19) 

We will now prove that a = 3/2, by examining the behavior of the generating function, 
which is numerically plotted on Fig. 1. Clearly, the radius of convergence of the series 
a(t) = J2i>i cat' (i-e. the location of its singularity) is t* = 1/R. At t = t*, a(t*) ~ J2i 1/'^" 
converges if a > 1, and its first derivative a'{t*) ~ J2i l/^"~^ diverges if a < 2. Then, Fig. 1 
indicates that 1 < a < 2. 

More generally, the behavior of a{t) around t* is given by a. Let us consider the singular 
part y{e) = a{t*) - a(t) with t = ^*e~^ For small e, y{e) ~ - e""]/z". For a > 2, the 
sum is dominated by small i and y{e) = 0(e), i.e. the first derivative a'{t*) is finite. For 
1 < a < 2, the sum is dominated hj i = 0(l/e) and y{e) = 0{e"~^). The case a = 3/2 
corresponds to y{e) = 0{y/e). 

Rewriting Eq. (|) 

ln[l + a(t)] = i[a(t)+t]+r(t) with r{t) = J2 ^Ht') + t% (20) 

^ k>2 ^'^ 

for t < t*, r{t) is a regular and increasing function. By defining g{x) = 21n(l + x) — x, 
t + 2r(t) = g{a{t)). The function g{x) has a maximum at x* = 1 with g{x*) = 2 ln2 — 1 and 
a parabolic shape (i.e. a finite and negative second derivative). As a{t) = g^^{t + 2r(t)), 
this gives a singularity for 

r + 2r(r) = 21n2- 1 and a(t*) = 1 (21) 
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with a square root shape a{t*)—a{t) — 0{^/t* — t). This proves that a = 3/2. Unfortunately, 
with this argument, we can not compute the vahic of t* = because of the complexity 
of r{t). Then R is sensitive to the details of the regular part of the recurrence : R is not 
"universal" . But a depends only of the kind of singularity : it is an "universal" exponent. 
As we have seen above, it is not surprising that it is equal to the a of the Catalan numbers. 

V. CONCLUSION 

We have seen that some concepts of statistical physics can be transposed to a problem 
of graph theory. The partition function of two-terminal series-parallel networks have an 
asymptotic behavior similar to those of self-avoiding walks or lattice animals. Its exponential 
growth depends on details on the model. On the other hand, a power-law correction, with 
an "universal" exponent can be computed exactly and have a fractional value —3/2. For the 
self-avoiding walks, this exponent is sensitive to the dimension and the kind of boundaries 
conditions and others independent exponents are related to geometrical properties. In our 
problem, it will be interesting to have a similar interpretation. For example, the number of 
terminal nodes (two in this article) could play the role of boundary conditions. By studying 
quantities like the average depth of elements (i.e. the number of decomposition needed to 
isolate this element), it is probable that new exponents would appear. 
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